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We consider the magnetophoresis problem within the Rubinstein - Duke model, i.e. a reptating 
polymer pulled by a constant field applied to a single repton at the edge of a chain. Extensive 
density matrix renormalization calculations are presented of the drift velocity and the profile of the 
chain for various strengths of the driving field and chain lengths. We show that the velocities and 
the average densities of stored length are well described by simple interpolating crossover formulae, 
derived under the assumption that the difference between the drift and curvilinear velocities vanishes 
for sufficiently long chains. The profiles, which describe the average shape of the reptating chain, 
also show interesting features as some non-monotonic behavior of the links densities for sufficiently 
strong pulling fields. We develop a description in which a distinction is made between links entering 
at the pulled head and at the unpulled tail. At weak fields the separation between the head zone 
and the tail zone meanders through the whole chain, while the probability of finding it close to the 
edges drops off. At strong fields the tail zone is confined to a small region close to the unpulled edge 
of the polymer. 

PACS numbers: 47.50.+d, 05.10.-a, 83.10.Ka 



I. INTRODUCTION 

The magnetophoresis problems is a member of the class 
of reptation problems in which a long polymer is driven 
through a gel. The reptative motion of the polymer can 
be succesfully modeled by a lattice version using sections 
of the polymer (the reptons) as the mobile units, which 
hop stochastically from cell to cell. The driving field 
is incorporated as a bias in the hopping rates, favoring 
the motion in the field direction. A simple and adequate 
model is the Rubistein - Duke model , which represents 
the polymer as a chain of N independently moving rep- 
tons, with the restriction that the integrity of the chain 
is preserved. The reptons trace out a connected string 
of cells in space, each cell containing at least one repton. 
The cells, which can be multiply occupied, carry the ex- 
tra reptons as units of stored length. In order to preserve 
the integrity of the chain, only those reptons which are 
located in cells with stored length can hop. 

Usually one considers the case of polyelectrolytes in 
which the reptons are uniformly charged. Thus the driv- 
ing field pulls equally on all reptons and the bias is the 
same for the hopping rates all along the chain. A practi- 
cal situation where this occurs is in DNA electrophoresis 
0- The DNA molecules, being acid, get charged in solu- 
tion and when they are placed in a gel subject to an exter- 
nal electric field they perform a biased reptative motion 
along the field direction. Electrophoresis is a technique 
of great importance in molecular biology and sequence 
analysis, as it allows to separate DNA strands according 
to their length 0. 

In nature the charge distribution is of course not al- 
ways uniform. The extreme alternative is the case where 
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FIG. 1: Examples of configurations of reptating polymers in 
the Rubistein - Duke model in the case of (a) Electrophore- 
sis and (b) Magnetophoresis. Black reptons perform a biased 
motion along the direction of the applied field while white 
reptons are unbiased. The configuration for a chain with TV 
reptons is given by a set of N — 1 integers (1/1,2/2 • • ■VN—i) 
measuring the distance of two neighboring reptons along 
the field direction (thus j/i = 0, ±1). For the two exam- 
ples shown the coordinates are (1, 0, 1, 1, 1, 0, 1, 1) for (a) and 
(1,1,1,0,-1,1,1,1,0,1,0) for (b). 



all reptons are neutral except one end repton, which is 
charged. A possible realization of such a situation is 
a magnetic bead, attached at one end of the polymer, 
which is driven by a magnetic field. Therefore this case 
can be referred to as the magnetophoresis (MP) problem 
It is the subject of this paper. As we will frequently 
compare our findings with the more common case of the 
uniformly charged polymer, we will, for briefness, refer to 
the latter as the electrophoresis (EP) problem and to the 
present one as the MP problem, although the distinction 
between the two is not electromagnetic, but only in the 
forces exerted on the reptons (see Fig. [Q. 
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The MP problem, within the framework of the 
Rubinstein-Duke model, has so far been studied by 
means of Monte Carlo simulations Q, and the calcula- 
tions were mostly restricted to the drift velocity as func- 
tion of the applied field and chain length. In this paper 
we analyze the MP problem by means of density-matrix 
renormalization-group (DMRG) techniques, which allow 
us to perform a detailed analysis of both global quantities 
as drift and curvilinear velocities and diffusion coefficient, 
but also on local average shapes of the polymer. 

The dynamics of the reptating chain is governed by the 
Master Equation which we put in the form 



(i) 



Here y stands for the set of links yi, ■ ■ ■ , Un-1, where yj 
measures the distance between the repton j and j + 1, 
along the direction of the applied field. In our lattice 
representation the yj can take the values ±1 and 0. The 
value yj = corresponds to the case that the reptons j 
and j + 1 occupy the same cell. Thus each zero is a unit 
of stored length. The non-zero values represent the cases 
where j + 1 occupies a cell "higher" (1) or "lower" (-1) 
than j. Higher and lower refer to a position in the direc- 
tion of the field, y represents a complete configuration 
of the chain (see Fig. 0. P(y, t) is the probability dis- 
tribution of the configuration at time t and the matrix 
H(y,y') contains the gain and loss transitions from y' 
to y. The bias in the hopping rate is contained in the 
matrix elements of H. Generally we have for the bias 
factor 



Bj = cxp{aqjE / ksT) 



(2) 



with qj the charge of repton j, E the driving field, a the 
distance between adjacent cells (measured along the field 
direction) and fcsP the standard combination of Boltz- 
mann's constant and the absolute temperature. In the 
MP problem all qj — except for j = N. We put 



B N = B = cxp(e/2) 



(3) 



and use e as the parameter for the driving field. The 
other parameter of the model is the number of reptons 
N. 

We have chosen this "hamiltonian form" of the Mas- 
ter Equation in order to stress the formal correspondence 
with a quantum mechanical model governed by a hamil- 
tonian matrix. The DMRG method exploits this analogy 
and indeed its success in one-dimensional quantum prob- 
lems carries over to reptation problems 0, 0, . We are 
interested in the stationary state of the probability dis- 
tribution. In the quantum language this corresponds to 
finding the right eigenvector of H belonging to the eigen- 
value zero. The adaptation of the DMRG method to the 
MP problem is straightforward and the data presented 
in this paper are obtained by the DMRG method. An 
important difference with the Master Equation is that in 



quantum mechanical problems the hamiltonian is hermi- 
tian, whereas in the reptation problem the matrix H is 
non-hermitian, due to the influence of the driving field. 
This restricts the applicability of the DMRG-method to 
moderately long chains and/or small driving fields. 

The physics of the EP and MP problems is qual- 
itatively different. As illustration consider the weak 
field case (e small), where the Nernst-Einstein relation 
v = FD relates the drift velocity v to the total applied 
force F and the diffusion coefficient D. The force F, 
equals Ne in the EP case, since one pulls at each rep- 
ton. It is well known that v scales as v ~ £/N , yielding 
for the diffusion the non-trivial result D ~ A~ 2 . In the 
MP problem D will be essentially the same, as hopping 
is limited by the availability of stored length. In both 
cases the motion can be considered as diffusion of stored 
length. Since F — e in the MP problem, we expect the 
drift velocity to scale as v ~ eN~ 2 , a feature which is 
born out by our calculations. 

In the Sec. |H]we discuss some moment equations de- 
rived from the Master Equation which are more helpful 
then in the EP problem in analyzing the drift- and curvi- 
linear velocity. They are expressed in terms of the prob- 
abilities rij that the link j has the value yj = k. The 
sum of the rfc adds up to 1 



(4) 



So it suffices to consider the two quantities n° and raj 
defined by 



(5) 



n°j can be called the local density of stored length, rrij is 
a measure for the local orientation and will be referred 
to as the profile of the chain. 

In the Sections IIIII and llVl we present our data for the 
velocities and the profiles. The analysis is most transpar- 
ent in the strong field limit where we can make an ansatz 
which almost perfectly represents the data. In section llVl 
we discuss the behavior of the profile for weak and strong 
pulling fields. 



II. MOMENTS OF THE MASTER EQUATION 

The DMRG method deals with the whole probability 
distribution -P(y). In the MP problem it is fruitful to 
consider moments of the Master Equation. One set of 
moments is obtained by multiplying Eq. Q with yj and 
then summing over all y. This leads to N — 1 relations 
which can be seen as an expression of the fact that the 
drift velocity v across all the N — 1 links of the chain 
must be the same on the average. An even more useful 
set of relations is obtained by multiplying the Master 
Equation by yj and summing over all y. The resulting 
N — 1 relations are an expression of the fact that the 
curvilinear velocity J is the same across all links. These 
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relations obtain the form 

J = - n§ 1< j < N (6) 

which can be viewed as the familiar law that the current 
J equals minus the gradient of the density of the stored 
length. In addition to Eq. JBJ one has two relations 01 
concerning the traffic in and out both ends of the chain. 
They read 

J = 1 - 3nJ, 

J = (n°r_i - + «_ a - n+^B' 1 (7) 

The expression for the drift velocity involves correla- 
tions between neighboring links 

v = ((l-y?_ 1 )y j -(l-y])y j - 1 ) (8) 

and it is therefore not as informative as Eq. JfJJl. As one 
sees Eq. (J7J involves only averages over the first (last) 
link. This holds also for the expressions for the drift 
velocity in terms of the averages of the first and last link. 

v = mi, 

« = (n° N _ 1 +n N _ l )B-{n Q N _ 1 +n + N _ 1 )B- 1 (9) 

These equations have been derived by Barkema and 
Schiitz |3§ using balance arguments. 

Equation is a powerful relation since it allows to 
express the density of stored length in terms of the curvi- 
linear velocity 

n° = n° - (j - 1) J (10) 

showing that the density profile is linear in the position 
of the cell. In particular IjlUfl implies a relation between 
the densities of the first and last cell 

n%-i = A -{N- 2)J (11) 

We have used the linearity of the density n° as a check 
of the numerical calculations. 

Counting the number of unknowns (v, J, n", n° N _ 1 , 
mi, rriN-i) and the number of equations (0), (0 and 
(|llfl we see that we have one more unknown than equa- 
tions. This situation is similar to the EP problem. There 
the expression for the curvilinear velocity does not obtain 
the simple form ©, due to the bias on the internal rep- 
tons. So one misses relation <|11|) . On the other hand 
J = for EP, due to the symmetry of the polymer on ex- 
changing head and tail. Thus in both cases the moment 
equations are not sufficient to determine the velocities. 
Higher moments do not lead to additional information 
since again higher order correlations appear. 

III. GLOBAL QUANTITIES 

We discuss first the behavior of global quantities as the 
drift and curvilinear velocities and the diffusion constant. 



-1 




log 10 (e) 

FIG. 2: Plot of logv vs. loge for various -/V. 

A. The Weak Field Limit 

In the weak field limit the polymer assumes mostly a 
random configuration and all the densities are close to 
1/3. The overal behavior of the drift velocity v as func- 
tion of e and N is given in Fig. [2J Note that the drift ve- 
locity becomes proportional to e for small e, as expected 
on the basis of the Nernst - Einstein relation discussed 
above. For stronger fields, the velocity saturates to a fi- 
nite values, as discussed in the next paragraph. A similar 
dependence of v on s and N was observed in the Monte 
Carlo study of Ref . |3J . The limiting behavior for small 
e and large N is thus: 

v(N) ~ sD(N) ~ ^, (12) 

with D(N) the zero field diffusion coefficient. The scaling 
behavior of D(N) as function of the length TV has been 
studied quite intensively 0, H, U El- Reptation the- 
ory predicts that D(N) ~ 1/-/V 2 , while conflicting results 
appeared on experimental measurements, for which both 
1/N 2 and, more recently, 1/A^ 2 3 jnj have been reported. 

A detailed study of the scaling of D(N), within the 
Rubinstein-Duke model by means of DMRG method was 
recently performed [^lil, for various end-point stretching 
rates. In that case the diffusion coefficient was calculated 
from the limiting value of the drift velocity for e — > 0, 
with the field acting on all reptons (the EP problem). 
Here we repeat the same analysis only for a single case 
(using a stretching rate d = 1 following the definition 
of d of Refs. 0,01). The advantage of calculating D(N) 
with a small field acting only on an end repton is that the 
DMRG procedure is much more stable in this case and 
one can compute longer chains. This is due to the fact 
that in the MP problem non-hermiticity is restricted only 
to the repton where the field is applied. As mentioned in 
the introduction non-hermiticity hampers the efficiency 
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FIG. 3: Plot of the effective exponent a(N) calculated for e = 
1CF 3 up to N = 51 from the decay of the diffusion coefficient 



D(N) and plotted as a function of 1/y/N. The fact that a(N) 
approaches linearly the limiting value 2 supports the scaling 
form for the diffusion coefficient given in Eq.(13). 



of the DMRG method. 

In order to calculate the diffusion coefficient from the 
Nernst-Einstein relation D — lim e ^o v /s m practice, we 
used a small field (e = 10~ 3 ) and checked explicitely 
that results do not change for smaller fields. The scaling 
behavior of the diffusion coefficient is expected to be: 



D(N)N 2 = A 



A' 



N 



(13) 



with A and A' some constants. The form of the sub- 
leading correction to D(N) has been debated for a while 
[ilj. Il5| and recent DMRG results suggest that it is of 
the type 1/y/N 0, supporting Eq ijTl)) . The coefficient 
was determined exactly HHHEg: A = 1/3. 

To analyze the scaling behavior of D(N) it is most con- 
venient to use the logarithmic derivative of the DMRG 
data: 



a(N) = - 



In [D(N)} - In [D(N + 2)] 
In N - ln( N + 2) 1 



(14) 



which is shown for N = 9, 11, . . . , 51 in Fig. |3J Plugging 
Eq. Ijl3(l in Eq. (|14(l one finds for the effective exponent 
a(N) — 2 + A' I (2Ay~N), a behavior which is accurately 
reproduced by our numerical data of Fig. [3] The present 
results corroborate previous claims |j| about the scaling 
form of D(N). 



B. The Strong Field Limit 

In the strong field limit the polymer assumes an ori- 
ented configuration, with the '+' links dominating at the 



FIG. 4: Difference between the drift v and curvilinear J ve- 
locities as function of the applied field and for various chain 
lengths. 



pulled end. At the other end we still have a substantial 
amount of links '0', since the polymer can only move by 
the diffusion of stored length from the tail to the head. 
Eliminating n\ from Eq. i|10|) with the use of Eq. JJJ we 
get 



»8r-i - j(l 



KJ), 



(15) 



where K = 3N — 5. In order that n^-i stays finite for 
N — > oo, the curvilinear velocity must vanish as 



J - K- 



N 



(16) 
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FIG. 5: Comparison between the DMRG data (symbols) and 
the crossover formula (solid lines) of Eq. 119H for the curvi- 
linear velocity J. 
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As one sees from <|15f) this limiting value is not sufficient 
to determine the limiting value of density which 
is sensitive to the corrections to (1161) . For the strong 
field limit it is useful to relate the drift velocity to the 
curvilinear velocity. With @ and JSJ) we get 



J + 2(n N _ 1 B-n%_ 1 B~ 1 ), 



(17) 



Now, if the polymer is fully stretched, v and J become 
the same. In Fig. 01 we have plotted the difference v — 
J as calculated by DMRG for various fields and chain 
lengths. We note that it is small for all values of e and 
N, in particular for strong fields and that this tendency 
is enforced for long polymers. That it also is small in the 
small field limit is a consequence of the fact that both 
quantities vanish in that limit. In order for the difference 
to vanish we must have 



l N-l 



B 2 n N _ 1 



(18) 



Now we may use this relation as the 6th relation, which 
enables us to make all the desired quantities explicit func- 
tions of e and N. We find for instance 



J{B) = v(B) 



B 4 - 2B 2 



1 



'■jv- 



-i(B) 



K{B A + B 2 + 1) + 3B 3 ' 

1 + B/K 
B 2 + 1 + W/K+ 1/B 2 ' 



(19) 
(20) 



This explicit field dependence is compared to the data in 
Fig. [5] and Fig. |B| The agreement is excellent in both 
cases. Note also that Eq. (|19|) is consistent with Eq. Ijl6|l 
and that it provides the proportionality coefficient. 

However, the crossover formulae (|19(1 - <|20(1 do not de- 
scribe the subtile dependencies in the limit of small fields 
B = cxp(e/2) — > 1. In this limit the drift velocity van- 
ishes as d ~ e, while one observes from Q19[l. that the 
curvilinear velocity vanishes as J ~ e 2 . For this reason, 
in the limit e — > 0, the crossover formula l|19f> predicts 
v ~ e 2 /(3N — 5), in disagreement with the correct scal- 
ing behavior of Eq. I|12fl • The strong field limit does not 
suffer from this problem. We note, for instance, that in 
the limit N — > 00 the saturation value of the velocity for 
B — > 00 is in agreement with the exact expression given 
in Ref. [|: v = l/(3N-5). 



IV. PROFILES 

Next we discuss some profiles, i.e. the local orientation 
mi = (yi) as function of the segment position along the 
chain. We consider N — 1 segments, thus iV-reptons with 
the charged one at head position N. Fig. shows a plot 
of rrii/e as function of the scaled variable (i — 1)/(N — 2) 
for chains of various lengths and at fixed field e = 0.001. 
This profile corresponds to the linear regime where the 
drift velocity scales as v ~ e. The notable feature is 
a symmetry between head and tail with respect to the 
center of the chain, although the magnetophoresis prob- 
lem is clearly asymmetric. This symmetry can be shown 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 
1/B 1/B 

FIG. 6: Comparison between the DMRG data (symbols) and 
the crossover formula (solid lines) of Eq. (I2( )t for n%_ 1 . 




0.0 0.2 0.4 0.6 0.8 1.0 
(i-l)/(N-2) 

FIG. 7: Average profiles {yi)/e for various lengths and for 
e = l(T 3 . 

fl8| to be strict in the weak field limit. It disappears for 
stronger values of the field as Fig. |S]shows, where profiles 
are plotted for e = 1 and various lengths N. 

In order to analyze the data further we also plot the 
individual probabilities for having a +, or — at the 
site i of the chain. For small values of e (not shown here) 
the curves are all near 1/3, with a slight excess of + links 
at the head and a depletion of — links. The densities of + 
and — links are monotonically increasing and decreasing 
functions of the position i. For intermediate fields e = 1, 
the densities are more interesting and in Fig. [5] we plot 
the values of n?, nf and n~ for N = 51. The linear 
behavior for n° is consistent with Eq. (|10|) . The curve 
for is monotonically increasing, but that for n~[ is not 
monotonically decreasing. 
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The qualitative behavior of the orientation profile can 
be understood by considering the "origin" of the non- 
zero links (yi = ±1) as has been introduced by Barkema 
and Newman 8]. In the MP problem more links are 
created at the pulled head than at the tail. They stream 
gradually down to the tail. We can keep track for every 
link yi — ±1, whether it is formed at the head or at the 
tail. After sufficient time the chain is divided into two 
zones: a head zone and a tail zone. They are separated by 
a small intermediate region with zeroes (we do not follow 
the origin of the zeroes). The zones remain separated 
because the yi = ±1, created at the head cannot cross 
the yi — ±1 created at the tail. The division between 
the two zones fluctuates in time and occasionally the tail 
zone disappears, while very rarely (particularly at large 
fields) the head zone vanishes. The larger the force on 
the head, the larger the asymmetry between the head and 
tail zones. We supplement these speculations by making 
an assumption on the ratios 



n{j) =ptU)/Pi U) 



(21) 



where pf(j) is the probability of finding a ± at i when 
the division is at j. We put 



n(J) = r h 

nti) = l 



for 
for 



j < i 
j > i 



(22) 
(23) 



Since in the tail zone there is no distinction between + 
and — we have set the ratio equal to 1. The idea, un- 
derlying this assumption, is that the + and — links are 
interlocked. So while moving in their zone their ratio can 
not change. 

At position i the average number of nonzero links (+ 
and — ) is equal to 1 — We introduce fi as the fraction 
of such nonzero links which are in the head zone. One 
can express the densities nf in terms of fi as 




0.4 0.6 
(i-l)/(N-2) 

FIG. 8: As in Fig. Qfor e = 1. 




0.0 0.2 0.4 0.6 0.8 1.0 
(i-l)/(N-2) 

FIG. 9: Plot of the average densities {nf}, {n°) and for 
e — 1 and N = 51. Inset: Blow up of the density n~ , showing 
a non-monotonic behavior with a linear increase as function 
of i by approaching the pulled edge. 



(/^ + (W^)[l-n°] 



(24) 



In both equations the terms proportional to fi are the 
contributions that i is in the head zone while the terms 
proportional to I — fi are the contributions from the case 
in which i is in the tail zone. We see that in the magne- 
tophoresis problem the situation simplifies, since we do 
not have to worry about the tail zone. It drops out when 
we consider the profile 



r h -l 
1 



(25) 



Thus the profile m t is, apart from the known factor 1 — 
n°, directly related to the fraction fi. The latter has 
a simpler interpretation. It starts out at i = 1 with a 
value nearly zero, since the head zone will only seldomly 
extend over the whole chain. It ends at i — N — 1 at a 
value very close to 1, since the tail zone will hardly ever 
extend over the whole zone. We can use this fact to tie 
the ratio rh to the end point values, discussed earlier, by 
considering (|25|l for i = N — 1 



>yN-i) = 



r h -l, 



r h + l' 

and solving for r^- It leads to 

1 - n° 

(Vi) = (VN-i)fj 1 ^ n o l 



(26) 



(27) 



JV-l 
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This form contains only /j as unknown. We can draw 
some conclusions from Eqs. (|25H and (|27l) for weak fields, 
as well as for long chains at stronger fields. 



A. Weak Fields 

For t^Owe may put 

r h = 1 + a h e 



(28) 



The function 1 — will approach the limit 2/3, so (J22 
becomes 



ah , 



(29) 



For the zero field limit of the profile we can take the zero 
field limit of fi. It has the property that head and tail 
become equivalent or 



fi = 1 - In 



l-i 



1 (30) 



The zero field limit of has been determined in Q by 
Monte Carlo simulations. We note that H30|) is consistent 
with the mentioned (18| symmetry in (yi). One should 
have ah — 2 in order that the profile becomes 2e/3 at 
the head, as is observed (see Fig. [J. This is perfectly in 
agreement with the value rh = B 4 ~ 1 + 2e for small e. 
Combining (|29|l and the first Eq. © we find that fx is 
a measure for the drift velocity v. According to (|12() fx 
should vanish as 1/3A 2 . This result has been derived in 
□ 

Another feature of Fig. [7| seems to be the collapse 
of the data on a single curve. Further data on longer 
chains show that the flattening-off at the ends of the 
chain shrinks with the size of the system and that the 
slope in the middle slowly decreases. This is another 
manifestation of the slow approach towards the asymp- 
totic behavior for large N. 

Note that if the division between the head and tail 
region were located with equal probability on all sites of 
the chain then one would have simply fi = i/N, which 
from Eq. (|29|l implies a linear profile. The profile of Fig. 
0is linear only at the center of the chain, while it strongly 
deviates from linearity close to the edges. This implies 
that the probability of finding the division between the 
head and tail regions is flat in the center of the chain and 
drops off at the chain edges. 



B. Long Chains and stronger Fields 

In this case the head zone will be dominant beyond a 
certain point (i.e. fi = 1 for i > zq ) in the chain, thus 



the division between the head and tail zones is expected 
to become localized close to the end of the chain which is 
not pulled. The curves in Fig. convincingly show this 
behavior. It is interesting to note that when /, — ► 1, Eq. 
(H3} becomes: 



+ n Oi Th 
nj = [l-Ti, 



= [1-n, ] 



r h + 1 
1 

r~ h + 1 



(31) 



It immediately implies that both nf and n~ are linearly 
increasing functions of i in the head zone, being a 
linearly decreasing function of i (see Eq. (!10[l L This 
explains the monotonic increase of n~ close to the pulled 
end shown in the inset of Fig. Note that using Eq. 
Q31JI one can estimate rn from the ratio of the slopes of 
nf and n~ in the head region. We find a ratio rh — B 4 
in agreement with our crossover formulae. 



V. DISCUSSION 

We have presented a series of numerical and analyti- 
cal results for the MP problem in the Rubinstein - Duke 
model, where a single reptating polymer is pulled by a 
constant driving field applied to one polymer end. We 
have shown that the numerical data for the drift- and 
curvilinear velocities can be quite well reproduced by sim- 
ple interpolating formulas following from the assumption 
that both velocities are equal in the limit of long chains. 
Indeed the measured differences are small which shows 
that the polymer is fairly stretched by the pulling force. 

We studied also local quantities, as the profiles, which 
provide information on the shape of the reptating chain. 
These are quite well understood using a representation in 
which the polymer is divided into a head and tail region, 
with different ratios of + and — links. At small fields the 
division between the two regions meanders trough the 
whole chain, and the probability of finding it close to the 
edges drops off. At strong fields the division gets localized 
close to the free end of the chain. Moreover some profiles 
show an unexpected non-monotonic behavior which has a 
simple interpretation in the interface picture. The precise 
shape of the profiles at weak fields close to the polymer 
edges, both at finite N and in the asymptotic limit N — > 
oo, will be discussed in details elsewhere [l8j . 
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